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Abstract. Our paper presents new methods for finding and testing of weak peri¬ 
odic variability of stellar objects developed for the purpose of detecting expected reg¬ 
ular light variations of magnetic chemically peculiar (mCP) candidates in the Large 
Magellanic Cloud. We introduce two new periodograms of the mCP star, BS Cir 
(HD 125630), appropriate for rotating spotted variables and compare the results with 
those obtained by the well-known Lomb-Scargle periodogram. The usage of peri¬ 
odograms and the testing of the significance of the found period candidates are demon¬ 
strated with two examples: the observed and simulated observations of the magnetic 
field of the mCP star CQ UMa (HD 119213) and the mCP candidate OGLE LMC136.7 
16501. Three newly developed tests of the periodic variability - the shuffling, bootstrap 
and subsidiary ones, are presented. We demonstrate that the found periodic variations 
known with Signal-to-Noise ratio larger than 6 can be approved as real. 


1. Introduction 

The chemical peculiar (CP) stars of the upper main sequence display abundances that 
deviate significantly from the standard (solar) abundance distribution. The existence of 
strong global magnetic field specifies a subset of this class, the magnetic chemically 
peculiar (mCP) stars. 

The periodic variability of mCP stars is explained in terms of the oblique rotator 
model, according to which, the period of the observed light, spectrum, and magnetic 
field variations is identical to the rotational period. The photometric changes are due to 
variations of global flux redistribution caused by the phase-dependent line blanketing 
and continuum opacity namely in the ultraviolet part of stellar spectra (Krticka et al. 
2007, 2012). The amplitude of the photometric variability is determined by the degree 
of non-uniformity of the surface brightness (spots), the used passband, and the line of 
sight. The observed light amplitudes are up to a few tenths of magnitudes, standardly 

In the Milky Way, we know of a statistically significant number of rotational pe¬ 
riods for mCP stars deduced from photometric and/or spectroscopic variability studies 
(Renson & Catalano 2001; Mikulasek et al. 2007). Nevertheless, also extragalactic 
mCP stars were found in the meanwhile. 

After the first photometric detection of classical chemically peculiar (CP) stars 
in the Large Magellanic Cloud (LMC) (Maitzen et al. 2001), a long term effort was 
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spent to increase the sample (Paunzen et al. 2006). Finally, we were able to verify our 
findings with spectroscopic observations (Paunzen et al. 2011). 

In this paper, we present the tools of the time series analysis developed for finding 
and testing of the expected periodic light variations of mCP candidates in the LMC 
(Paunzen et al. 2013). The list of targets (Paunzen et al. 2006) was compared with the 
OGLE database (Udalski et al. 2008) for corresponding measurements. In total, four¬ 
teen common objects were found and the V and I light curves analysed. The description 
of methods arc also published in the textbook by Mikulasek & Zejda (2013). 


2, The periodograms 

The basic tool to search for periodic variations of a measured value (intensity, mag¬ 
nitude, and radial velocity) arc so called periodograms. These plot some quantities 
expressing the significance of the phase sorting of searched data according to an indi¬ 
vidual angular (circular) frequency to = 2n f = 2n/P, where / is a frequency, / = P ] , 
P is a period. The extremes of periodograms then indicate the best arrangement of ob¬ 
served data versus a period and hence the possible periodicity of a signal. 

The basic linear regression models of periodic variations arc harmonic polynomi¬ 
als of the g-order: 

g 

F(t, to) = ^ /?i/(m) cos(mj t) + fajioj) sin(wj t), (1) 

j= i 

where Fit, to) is the model of detrended measured quantities {}’,(?,)} corrected for their 
mean, fi\j(to), (hjito) are 2g harmonic coefficients. The harmonic coefficients for the 
best fit of model function F{t,to): by(to), byioo) for the fixed to can be determined by 
the standard least square method technique allowing to count with uneven uncertainties 
of individual measurements {cr,}. 


2.1. Periodograms with modulated amplitude 

The simplest way how to construct LSM spectral periodogram is to plot scalar value 
X~ito) versus to or / = P _1 , where 
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Now we can find and discuss the frequencies for which the value ;y 2 (o») reach their min¬ 
ima. This method is fairly general because it can be applied to any kind of time series 
(magnitudes, intensities, spectral line equivalent widths, or radial velocities). Never¬ 
theless for data of the same type (magnitudes, intensities) we recommend to use the 
following modification with some ’value added’. 

The first sum of equation (2) where the first sum on the right is a constant that not 
depends on the frequency, while the second is the weighted sum of the squares of the 
model prediction given by the function F(to, t). Therefore, instead of the minimum of 
the x 2 (° J ) we can And the maximum of the modulated amplitude A m 
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which is nearly equal to the effective amplitude A e ff(®) of a periodic function (see in 
Mikulasek et al. 2007). 

For the first estimate of the variation period it is sufficient to use the simplest 
possible option: g = 1 which gives also the most explicit results. Then 
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Figure 1. All three periodograms of the moderately cool chemically peculiar star 
BS Cir (HD 124224; for more information see e.g. in Mikulasek et al. 2013, 2015) 
undoubtedly pinpoint the only dominant period peak at P — 2 d 2042. Periodograms 
are displayed in arbitrary units normalized to the height of the period peak. The 
results of particular periodograms are comparable, nevertheless both of the new tools 
of periodograms ((b) and (c); Sec. 2.2 and Sec. 2.1) are a slightly better than the 
results of the notorious Lomb-Scargle periodogram ((a). Sec. 2.3). 


2.2, Robust S/N periodogram 

The second LSM type of periodograms uses as a criterion of the significance of indi¬ 
vidual peaks a robust “signal-to-noise” S/N criterion which is defined as: 


S/N(oj) = 


QM 

SQicoY 



(5) 
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Figure 2. The periodogram of the square of the modulated amplitude, constructed 
on the basis of 4766 individual photometric observations in the spectral interval 505 
to 560 nm, indicates the main frequency /o = P ~ 1 = 0.4538 d (rotation period) and 
subsidiary peaks at/o/2, 2/o, 3/o, 4 /q. 


5Q{co) is an estimate of the uncertainty of the quantity Q(u>) for a particular angular 
frequency. The model function F(u>, t, bn, bi\, • • •) is assumed in the form of the har¬ 
monic polynomial of the g-th order (see Eq. 1). The detailed description of both LSM 
novel periodogram criteria can be found in the textbook by Mikulasek & Zejda (2013). 

We tested the properties of the S /N(u>) criterion on thousands samples with sine 
(g = 1) signals scattered by randomly distributed noise. We found that if there is no pe¬ 
riodic signal in such data, the median of the maximum S /N{co) value in a periodogram 
is 4.52, in 95% of cases we find a S/N value between 4.2 and 5.4. Consequently, the 
occurrence of peaks definitely higher than 6 indicates possible periodic variations. 

The periodogram S/N{u>) is very sensitive (see Fig. lb) and suppresses the false 
periods which results from usually bad time distribution of astrophysical observational 
data, very well. 


2.3. The classical Lomb-Scargle periodogram 

During the treatment of OGLE-III time series (Paunzen et al. 2013), we concluded that 
both types of periodograms correlate very well with other time-proven periodograms as 
e.g. the Lomb-Scargle (see e.g. Press & Rybicki 1989) one. So we are able to consider 
them as generally interchangeable (see also Fig. 1). 

The Lomb-Scargle method assumes that the found changes have sine/cosine type 
phase curves. The method uses, as the measure of individual periods significance, the 
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quantity Qis(oj): 
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Z/Li cos 2 [w (ti - r)] + £"=1 sin 2 [w (f,- - r)] 


The definition of the significance indicator <2ls (see Eq. 6) can be easily modified to 
the more general case of data with uneven uncertainties. 


3. Examples of strong and weak variability 

3.1. Light variation of BS Cir - a well observed object with strong light variations 

BS Cir (HD 125630) is a moderately cool mCP star with the rotation period P = 2 d 24. 
There arc many observations of the double-wave light curves available (for details see 
in Mikulasek et al. 2013, 2015). Although light curves taken in various colours differ, 
those ones obtained in the spectral region 500-560 nm can be considered as almost 
identical. 

The periodogram 1 presented in Fig. 2 shows the spectrum of the square of the 
modulated amplitude A 2 //), (g = 2, see Eqs. 1,3) dominated by the basic period with 
frequency /o = 0.45363 d _1 , and subordinate peaks at fo/2, 2/o, 3/o,4/o- This spec¬ 
trum results from the fact that the light curves in various filters arc not exactly the 
same. The mean modulated amplitude of the maximum peak correspond to 0. 11 mag, 
the mean weighted S /N of used 4766 individual photometric observations is about 20. 

We obtain such a periodogram only in very rare cases when we study an object 
with the pronounced variability documented by thousands of observations not suffering 
from bad sampling. Nevertheless, frequently we have to analyze data of much worse 
quality and time sampling. 

3.2. Magnetic field of CQ UMa - strong variability, but few observations 

CQ UMa (HR 5153 = HD 119213), is also a moderately cool mCP star, similar to BS 
Cir, with a strong variation in the v color and c\ index. The period of variations of all 
kind is constant, P = 2 d 449967 Mikulasek et al. (2011). 

In 1983-4, we obtained nine measurements of the effective magnetic field (B e ff) 
using the hydrogen magnetograph of the 6-m SAO reflector (see full dots in Fig. 3). 
The expected synchronous periodic variations in the light and magnetic field (B e ff) were 
confirmed (Mikulasek et al. 1984). 

However, arc these nine magnetic observations of CQ UMa sufficient for the pe¬ 
riod determination? Definitely not, the period spectrum is overcrowded by plenty of 
aliases (Fig. 4). Several of them arc so much dominant that they exceed even the peak 
of the real period. 


’The periodograms of BS Cir were based on 4766 individual Hp, y and V measurements taken from Vogt 
& Faundez (1979); Manfroid & Renson (1983); Catalano & Leone (1993); ESA (1998); Pojmariski et al. 
(2001); Mas-Hesse et al. (2004); Mikulasek et al. (2013, 2015). 
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Figure 3. The phase diagram of real and artificial (simulated) magnetic field ob¬ 
servations CQ UMa. The periodogram of based on 9 measurements of its magnetic 
field obtained during 2 years. Real rotational period is denoted by the arrow. 


The aliases can be suppressed by a continuation observations. We can simulate it 
by adding of 27 ‘new’ (artificial) magnetic measurements obtained during the follow¬ 
ing six years (see open circles in Fig. 3). The new periodogram (Fig. 5) undoubtedly 
indicates the only dominant period at the true placement 2 . The phase diagram plotted 
with this period displays simple sinusoidal variations with a S/N ~ 8. 

If we remove the sinusoidal signal from our magnetic measurements, we can 
search for other possible relic periods using the same periodogram technique - see the 
bottom of Fig. 4. It seems that there is no periodic signal in the remaining data meaning 
that the found period is unique. The remaining data correspond to pure noise which is 
not very prominent. The characteristics of the periodogram is determined mainly by 
aliases originated from the very low number of measurements (36) and the short time 
basis of observations (8 years). 

The aliases can be lowered if we use, besides not very accurate magnetic mea¬ 
surements of Z? e ff, more precise photometric data (S /N usually up to 20). Then we arc 
able to achieve the pure periodogram similar to the periodogram of BS Cir displayed in 
Fig. 2. 


3.3. Is the mCP candidate OGLE LMC136.7 16501 periodically variable or not? 

Several quite different problems were encountered when we analyzed the light curves 
of mCP candidates in the Large Magellanic Cloud. The number of individual OGLE 


-The other dominant peaks in the periodogram are aliases conjugated with the period of one sidereal day 
through the Tanner relation - see Eq. 7 
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Figure 4. The periodogram of CQ UMa based on 9 measurements of its magnetic 
field obtained during 2 years. Real rotational period is denoted by the arrow. 


observations in the V and I colours are ~ 300 - 400 and the duration of observations 
(8 years) were satisfactory. The effect of aliases in such data should be insignificant. 
Unfortunately, the periodograms suffer here from relatively large scatter due to obser¬ 
vational noise. The ratio S /N of an individual observation was very low because of the 
weakness of the signal. 

We shall now discuss the period analysis of a typical representative of the mCP 
candidates in LMC - the star OGLE LMC136.7 16501 (04 50 46.10 -69 59 16.7, 
2000.0), denoted in Paunzen et al. (2013) as a star No. 4, with a mean magnitude 
in V = 19.16 mag, the absolute magnitude My = 0.41 mag. The periodogram depicted 
in Fig. 6 is quite dissimilar to the periodograms of mCP stars discussed above (compare 
with Fig. 2 and Fig. 5). We do not see here any prominent period peak. The characteris¬ 
tics of the period spectrum is mainly determined by the stochastic scatter that competes 
with several possible periods of the the principal aliases conjugated with the basic sam¬ 
pling period of one sidereal day. The frequencies of aliases obey to notorious Tanner 
(1948) relation: 

f a \(k) = |/m + kf s \, where k = 2, -1,1,2,...), (7) 

where f d \(k) is the frequency of the alias of the k-th order, /m is the frequency of the 
maximum period peak (/m = f a i(0)), and f s is the sampling frequency, here it is the 
frequency of a sidereal day / s = 1.00274 d -1 . 

The most significant period in the periodogram (Fig. 6) is 4/373, with amplitudes 
of the light curves A m = 0.044 mag, both in V and I. The ratio of S/N = 5.1 for this 
peak is very low. The light curves in V and I arc depicted in Fig. 7. Nevertheless, all 
predicted aliases: / a i(-2) = 1.777, / a i(-l) = 0.774, / a i(l) = 1.231, /ai(2) = 2.234 are 
indicated. 
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Figure 5. The periodogram of CQ UMa based on 27 measurements of its magnetic 
field obtained during 6 years. Real rotational period is denoted by the arrow. The 
dark periodogram corresponds to the period analysis if we subtract the variations 
with the basic period P = 2 d 449967. The low amplitude of it and the absence of 
any other pronounced period peaks mean that the periodogram of the star is formed 
namely by aliases of the basic (true) period. 


Does it mean that the period 4*?373 is real? Unfortunately not. The occurrence of 
aliases at their ‘right positions’ given by the Tanner relation Eq. 7 is only the confirma¬ 
tion of the distinctiveness of expected sampling period or periods in observational data. 
Because the OGLE data were obtained at a ground based station, the pronounced one 
day sampling period can be easily explained. The observed ordering of observation into 
a phase curve with the period 4 d 373 can be a mere coincidence. 

That is why we shall solve the crucial question: Are the found periodic variations 
real? Some partial answers may give us several simple tests of variability of general 
use we have developed (Paunzen et al. 2013). 

3.3.1. Shuffling method 

The heuristic shuffling method is able to test the ‘pessimistic’ hypothesis that the dis¬ 
tribution of the data is only random. We can analyze the data of our object displaying 
weak periodical variations again by the same way as it was described above, only we 
randomly shuffled all the individual observations (magnitudes and their uncertainties), 
the times of the observations remained the same. We are convinced that this derogation 
will destroy any periodic signal and randomize the data. For the shuffled version of data 
we then find the maximum modulated amplitude A m and the maximum of S /N. Then 
we repeat the same procedure many times and compare the results with A m and S /N for 
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Figure 6. The periodogram of LMC 136.7 16501 displays many period peaks. 
The period spectrum corresponds to nearly pure scatter. 



the original unshuffled data set. If the results are nearly the same, we may conclude that 
the periodic variability of the object (if any) is undetectable in the investigated data. 

We found for the mCP candidate OGLE LMC136.7 16501, the median of ampli¬ 
tudes in period peaks of shuffled data A m = 0.040 mag and a ratio S/N = 4.6, what is 
very close to the values found above A m = 0.044 mag, both in V and I and S/N - 5.1. 
Consequently, the LMC136.7 16501 photometric data are very probably random. 

The same results we obtained also for the other 11 mCP candidates, only for stars 
No 12 and No 14 (according to the list of LMC mCP cadidates in Paunzen et al. 2013) 
the periods seem significant. We deduce a S /N > 6 for significance. 

The following tests will help to quantify the significance of the periods. 


3.3.2. Bootstrap test 

The technique of bootstrap (Hall 1992) has proved to be very useful for testing the 
statistical significance and therefore the reality of found periods. It helped us to quantify 
this reality as a probability that the periodogram of a randomly created bootstrap subset 
has its dominant peak at the same frequency as the standard periodogram. We tested it 
with one hundred of bootstrap subsets for each star of our sample. We consider a period 
as statistically significant, if the maximum peak occurs at one of the aliased frequency 
because during the bootstrap choice aliases often exceed the basic peak. 
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Figure 7. I and V light curves of OGLE LMC136.7 16501 plotted versus the ‘best’ 
period P = 4 < ?373. 


The results for found periods of all 14 mCP candidates were dismal: the bootstrap 
significance never exceeds 50%. Especially, the bootstrap significance of the period 
4 d 37 in the case of the star 4 was only 26%! The periodicity of stars No. 12 and 14, 
mentioned in Sec. 3.3.1, is disputable. 

3.3.3. Subsidiary test 

For data, including observations in two or more spectral regions - their individual pe- 
riodograms should indicate nearly the same peak period. It can be tested, e.g. by the 
periodogram of the product of individual A m or S /N values. 


4. Instead of a conclusion 

All mentioned newly developed methods arc open source. For details see in Paunzen et 
al. (2013) or the textbook of Mikulasek & Zejda (2013). 
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